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Several lattice calculations of a combination of the low energy constants of the chiral Lagrangian, 2as — 0:5, 
are presented. This combination is critical for the preclusion of a massless up quark. The result found is 
2a8 — as = 0.115 ± 0.051"*"' ± 0.25'^-^'^*, which is well outside of the range allowed by a massless up quark. 



1. Introduction 

The symmetries of QCD allow for a CP break- 
ing term in the QCD Lagrangian. However, mea- 
surements of the neutron dipole moment have 
demonstrated that the coefRcient of this term is 
very small, if not zero. The unnatural smallness 
of this coefficient is known as the strong CP prob- 
lem. A massless up quark has long been proposed 
as a potential elegant solution to the problem. 
Without the up quark mass term, one is free to 
remove the CP violating term through a field re- 
definition. 

At tree level Chiral Perturbation Theory 
(ChPT) appears to answer the question of a mass- 
less up quark. The quark mass ratios dictate 
the form of the Lagrangian's mass matrix, which 
at lowest order fully determines the light meson 
mass ratios. 

However, the NLO terms in the chiral La- 
grangian contribute to these ratios as well. If this 
contribution is strong enough, it would mimic the 
effects of a massive up quark and allow a mass- 
less up quark to be consistent with experimental 
results. This is known as the Kaplan-Monohar 
ambiguity fltl. 

Distinguishing between a light and a massless 
up quark requires knowledge of the coefficients 
of the NLO terms in the chiral Lagrangian, the 
Gasser-Leutwyler (GL) coefficients. Specifically, 
it is the combination of constants 2a8 — '350 which 
corrects the relevant ratio. If this combination 
falls within a certain range, —3.3 < 2as, — as < 



1 s 




*The low energy constants on are related to the corre- 
sponding constants of the chiral Lagrangian Li by the re- 
lation ai = 8{4:n)^Li. 



Figure 1. 16^ x 32, /3 = 5.3, mg = 0.01, hyp; 2a8 - 
Q5 = 0.115(51). The x's were excluded from the fit. 

— 1.5, current experimental results can not rule 
out a zero up quark mass. 

The full set of GL coefficients, including the 
combination of interest, can not be determined 
by experiment. Various phenomenological argu- 
ments have been used in the past to assemble a 
standard set of values for the coefhcients. How- 
ever, these coefficients are determined by the low 
energy non-perturbative behavior of QCD, and 
thus the lattice offers the best opportunity for 
calculating them directly. 

2. Partially Quenched Chiral Perturbation 
Theory (PQChPT) 

PQChPT is the tool with which one can calcu- 
late the GL coefficients on the lattice. PQChPT 
is distinct from standard ChPT in that it is 
constructed from the symmetry of a graded 
group. This graded group follows from the pre- 
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Table 1 

Simulation details. 



L 


T 


/3 


ms 


Nf 


a"-' (MeV)f 


hyp* 


2as — as 


16 


32 


5.3 


0.01 


3 


1347.9(31) 


V 


0.115(51) 


16 


32 


5.3 


0.01 


3 


1260(28) 




-0.071(14) 


8 


32 


5.3 


0.01 


3 


1366(49) 




0.47(12) 


8 


32 


5.115 


0.015 


3 


730(130) 




0.310(19) 


8 


32 


5.1235 


0.02 


3 


703(97) 




0.314(17) 


8 


32 


5.132 


0.025 


3 


780(120) 




0.375(13) 


8 


32 


5.151 


0.035 


3 


800(130) 




0.437(18) 



sumed quark content of Partially Quenched QCD 
(PQQCD), separate valence and sea quark flavors 
in addition to ghost quark flavors, which in per- 
turbation theory cancel loop corrections due to 
valence quarks. 

The Lagrangian of PQChPT up to 0{p^) fol- 
lows, with only relevant NLO terms shown. 

C ^ ^ sTr [^^,U^^'U^] + ^— sTr [xU^ + Ux] 

+ Uslr[d^Ud^'U'^] slr[xU^ + Ux\ 
+ L5 sTr [5^C/9^C/t {xU^ + Ux)] 

+ L^slr[xU^ + Ux]' 

+ Ls sTr [xU^xU^ + UxUx] + ■■■ (1) 

U = exp (2^$//) (2) 
X = 2/xdiag({ms,my}) , (3) 

where $ contains the pseudo-Goldstone "mesons" 
of the broken SU{Nf + Nv\Nv)l » SU{Nf + 
Nv\Nv)r symmetry and U is an element of that 
group. In our calculations three degenerate sea 
quarks were used, Nf = 3, while the number of 
valence quarks Ny cancels in all expressions, af- 
fecting only the counting of external states. The 
constants /, /.t, and Li are unknown, determined 
by the low energy dynamics of PQQCD. 

It was recently realized that, because the va- 
lence and sea quark mass dependence of the 
PQChPT Lagrangian is explicit and because full 
QCD is within the parameter space of PQQCD 
(my = ^,5), the values obtained for the GL co- 
efficients in a PQQCD calculation are the exact 
values for the coefhcients in fuU QCD Fur- 
thermore, the independent variation of valence 
and sea quark masses allows additional lever arms 



in the determination the these coefficients. Be- 
cause the Nf dependence of the Lagrangian is not 
explicit, the GL coefficients are functions of Nf. 
Thus, it is important to use a physical number of 
sea quarks, as we have, when extracting physical 
results. 

3. Predicted Forms 

PQChPT predicts forms for the dependence of 
the pseudoscalar mass and decay constant on the 
valence quark mass, here assuming degenerate sea 
quarks and degenerate valence quarks, and cut- 
ting off loops at — 47r/. 

ml = {4:Trf)'^zinv ll + zmv (ias - as + ^ 



— (2my - ms) In zmy 
Nf 



(4) 



zmy 



zNf z I 

+ [mv + ms)l-D--[mv + ms) I (5) 



z = 



2^1 



(6) 



These forms differ slightly from those in as 
the NLO dependence in the sea quark mass has 
been absorbed into fi and /. This is allowed as the 
error due to this change manifests when z appears 
in the NLO terms, pushing the discrepancy up 
to NNLO. Accounting for these absorbed terms 
would require a systematic study at several sea 
quark masses. 

The forms above are derived assuming degen- 
erate light mesons. However, our use of staggered 
fermions on rather coarse lattices is likely to gen- 
erate significant flavor symmetry breaking, and 
thus a splitting of the light meson masses. Al- 
though there is no simple quantitative method to 
account for this split, one effect might be a reduc- 
tion in the importance of meson loops, and thus 
a weakening of the chiral log terms. In order to 
study the systematic errors due to ffavor symme- 
try breaking, we applied hypercubic blocking Q 
to our primary ensemble. 



t Lattice spacing determined via rg. 
Denotes a hypercubic blocked ensemble. 
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Figure 2. 8^ x 32, several values of (3 and mg with 
matched lattice spacing. 

To determine the value of 2as — 0.5, the local 
pseudoscalar correlator was calculated using sev- 
eral valence quark masses. The correlators were 
fit to exponentials, while simultaneously the me- 
son mass and decay constants were fit to the pre- 
dicted forms. For the smaller 8'^ x 32 ensembles, 
a constant term was added to the meson mass 
form, consistent with finite volume. The results 
of the fit are the values /i, /, a^, and 2as — a^- 

Figure |l| and || display the quantity 



Rm 



Tnl{ms)mv 



(7) 



suggested by [Q. This quantity accentuates the 
NLO terms in mj, and thus is useful for plotting 
purposes. It should be noted, however, that the 
full forms of and were used when fitting. 
When calculating the ratio, we did not use the 
simplification shown in 0| , but rather used a full 
numerator and denominator. The value ttik in 
the figures is the valence quark mass at which 
the meson mass equals the physical kaon mass. 

4. Results 

2a8 — as was found to vary significantly among 
the ensembles. Changes in volume altered the re- 
sult, suggesting that yet larger volumes may be 
required. Additionally, the application of hyper- 
cubic blocking shifted the result, suggesting that 
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Figure 3. Compiled results. The UKQCD Collab- 
oration data point is taken from Q, showing their 
statistical error only. The rightmost point shows the 
range allowed by a massless up quark. 

flavor symmetry breaking is indeed affecting the 
value. Despite systematic shifts in the calculated 
value between ensembles, it remained well outside 
the allowed range for a zero up quark mass. 
The quoted result of 2as — — 0.115 ± 

0. 051.*^* ± 0.25''^*'* comes from our hypercubic 
blocked 16'^ x 32 ensemble, with generous sys- 
tematic errors due to the significant fluctuations 
in the result between ensembles. This falls out- 
side the range allowed by a zero up quark mass, 
-3.3 < 2ag, ~a5 < -1.5. 
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